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We study the dipole collective oscillations in bose-fermi mixtures in a dynamical time- 
dependent approach, which is formulated with the time-dependent Gross-Pitaevskii equation 
and the Vlasov equation. While the bose gas oscillates monotonously, the fermion oscillation 
shows a beat and damping. We find big differences in behaviors of fermion oscillation between 
the time-dependent approach and usual approaches such as the sum-rule approach. When 
the amplitude is not minimal, the dipole oscillation of the fermi gas cannot be described 
with a simple center-of-mass motion. 
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I. INTRODUCTION 



Over the last several years, there have been significant progresses in the production of ultracold 
gases, which realize the Bose-Einstein condensates (BEC) |],y,y,B^two boson mixtures B, 
degenerate atomic Fermi gases 6], and Bose-Fermi (BF) mixing gases p, yjl^- In particular the 
BF mixtures attract physical interest as a typical example in which particles obeying different 
statistics are intermingled. Using this system we have a very big opportunity to get new various 
knowledge about many body systems because we can make a large variety of combinations of 



atomic species and control the atomic interactions using the Feshbach resonance 



excitations 
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studies of the BF mixtures have been done for static properties 111 
phase diagram and phase separation 17, for stability 
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Above all phenomena the spectrum of the collective excitations is an important diagnostic signal 
for these systems. Such oscillations are common to a variety of many-particle systems and are 
often sensitive to the interaction and the structure of the ground state and the excited states. 
Theoretically collective motions are usually studied with the random phase approximation (RPA) 



28( 1 or its approximate methods such as the sum-rule [2S] and the scaling 



approaches. The first author (T.M.) and his collaborators [3]| studied monopole oscillations 
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of BF mixtures by calculating the time evolution of these oscillations with the time-dependent 
Gross-Pitaevskii (TDGP) equation and the Vlasov equation. This dynamical approach showed 



different behaviors of the oscillation from RPA 



281 ]. such as rapid damping at zero temperature. 



The RPA can treat only states with one-particle one-hole excitation and only describe minimal 
vibrations around ground states; for example it is shown in Ref. 28|] that the radial variation is 



only about 0.05 % in RPA, while the amplitudes of actual experiments can be as large as 10 % 
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For the single boson syst 



collective motions well 



em the RPA 35| explained experimental results on frequencies of 



341 ]. In this system all bosons occupy one single particle state at zero 



temperature, and their collective motions are simple; damping does not appear at very low 
temperature in experiments |3(1] and in a theoretical work with the time-dependent density 
matrix theory [s^]- 

On the other hand fermions occupy many single particle states even at zero temperature. When 
the amplitude is about 10 % of RMSR, the monopole states have as large excited energy as 
ten to one-hundred times of the one-particle and one-hole excitation energy. Their collective 
oscillations with a large amplitude are multiparticle and multihole states and include various 
mode with different frequencies. Indeed an oscillation of population difference in two-component 
fermi gas has shown damping due to a multimode dephasing 83]. In the BF mixtures, especially. 



fermions occupy a larger region than condensed bosons, and have different potentials between 
inside and outside of the boson occupation region. Then their motions are not harmonic, and 
then the fermion oscillation makes damping [31]. In order to study the collective oscillation in 
BF mixtures, hence, we need to calculate time evolution of the system using a time dependent 
dynamical approach. 

Thus the BF mixing gases show new dynamical properties different from those in other finite 
many-body system such as nuclei. In order to have more information on the dynamical proper- 
ties, we need to investigate other kinds of multipole motions. In this paper we study the dipole 
oscillations in bose-fermi mixtures as the next step by solving the time evolution of the con- 
densed boson wave function and the fermion phase-space distribution function with the TDGP 
equations and the Vlasov equations, respectively. 

In the next section we explain our transport model to calculate the time evolution of the system. 
In Sec. Ill we show the calculational results for the dipole oscillation in the BF mixture, and 
discuss their properties. Then we summarize our work in Sec. IV. 
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II. TIME EVOLUTION EQUATIONS 

Here we briefly explain our approach. In this work we consider a dilute boson and one- 
component-fermion coexistent gases at zero temperature with the axial symmetry with respect 
to the z-axis. We assume only zero-range interaction between atoms, and there is no fermion- 
fermion interaction in the system. The hamiltonian is written as 



H = 



I 



2Mb 

re .„oT. . 1 



Mr 



(1) 



where (j) and ■0 are boson and fermion fields, respectively, Mb and Mp are the boson and 
fermion masses, Mj. = MbMp/{Mb + Mp) is the reduced boson-fermion mass, CIb and Clp are 
the transverse trapped frequencies of the boson and the fermion, ubb and ubf are the s-wave 
scattering lengths between two bosons and between boson and fermion, respectively. In addition 
the positional coordinate is described as g = (qt^Ql), and kl is the ratio of the longitudinal 
trapped frequency to the transverse trapped frequency. 

In this formulation we can change all variables to dimensionless ones without losing generality 
as follows. We normalized the spatial coordinate q to the dimensionless coordinates r as r = 
{rT,z) = q/Rb with Rb = {Ti/Mb^bY^'^ ■ According to this normalization the boson and 

— 1/3 ~ ~ 

fcriniou fields are also scaled as ^(V') = Rb <A(V')- By dividing the above hamiltonian H by 

KIb, then, we can define the dimensionless hamiltonian H = H/hO,B as 



H 



2mf ' ' "^^^ ' 2 
+hBF<t>\r)<t>{r)il^\r)il^{r)\ , (2) 

where nif = Mf/Mb, cof = ^f/^b, and gBB = '&'nhaBBR^ and Ubf = 4:TrhmfaBF{^ + 
mfr^Rs'- 

In this work we consider the zero-temperature system, including N^, bosons and Nf fermions, so 
that the total wave function is written as 

'^ir) = {llMr^)}^f[i'n], (3) 

i=l 

where (/)c is a wave function of the condensed boson and is a Slater determinant of fermions 
with single particle wave functions, V'n' 
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The time evolution of the wave functions are obtained from the variational condition that 

dr < $(r)|{i— -F}|$(r) >= 0. (4) 
From this condition we derive coupled equations of the TDGP and TDHF equations as follows. 

i^Mr,T) = |-^V2 + C/B(r)| </.,(r,T), (5) 



ij-Mr,r) = -i^^i- + Upir) ^r^r) (6) 

OT I ZlTLf J 

with 

UB{r) = ]^{r% + k\z^) + gBBPB{r) + hBFPF{r), (7) 
Upir) = ^mfUjj{rT + Klz^) + hBFPB{r), (8) 
where ps and pp are boson and fermion densities which are given as 

PB{r) = m\Mr)\\ (9) 

occ 

PF{r) = j;iVn(r)p. (10) 

n 

The number of fermion states are usually too large to solve the above TDHF equations directly 
so instead one uses the semi-classical approach. In the semi-classical limit {h — > 0) the TDHF 



equation is equivalent to the following Vlasov equation 



^fir,p;T) = ||- + —Vr - [VrUFir)][vAfir,p;T) = 0, (11) 
ar lOT nif J 

where f(r,p;T) is the fermion phase-space distribution function defined as 

f{r,p,T)= [ d\<m^{r + -u,T)ij^{r--u,T)\<l>>e-'P^. (12) 
J 2 2 



As an actual numerical method we introduce the test particle method 



40l | to solve the Vlasov 



equation (jlip and describe the fermion phase-space distribution function as 

f{r,p,T)=^-^ 5{r-r,{T)}5{p-p,{T)}. (13) 
where Nt is the number of test-particles per fermion. 

By substituting Eq. (|13|) into Eq. (|ll|) . we can obtain the following equations of motion for test- 
particles as 

|:ft(T) = -V,,(7F(r), (15) 

Thus we obtain the time evolutions of the condensed boson wave function and the fermion phase 
space distribution function by solving eqs.([S]), (HH) and (fT3|) . 
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III. RESULTS AND DISCUSSIONS 

In this section we show the results of our calculations on the dipole oscillations in the BF 
mixture. We deal with the system ^''Rb -^''K, where the number of the bosons C^^Rb) and the 
fermions (^'^K) are A*";, = 10000 and Nf = 1000, respectively, and the fermion mass normalized 
by boson mass is ruf = 40/87 ~ 0.46. We assume the spherical trap {kl = 1) and take the 
fermion trapped frequency to be = 1/ J ruf ~ 1.48. The boson-boson interaction parameter 



qbb is fixed to be qbb = 1-34 x lO"'^ [4]|, which corresponds to qbb = 4.22(nm), while the 
BF interaction parameter Kbf is varied. In the numerical calculation we take the number of 
the test-particles per fermion to be = 100 and solve the time evolutions by using the second 
order predictor-corrector method. 

A. Ground State 



The sum-rule approach [29(] showed that the dipole oscillation frequency is sensitive to the density 
distributions. In this subsection, then, we see the boson and fermion density distribution from 
the point of view of the boson-fermion coupling dependence. 

In the ground state the the wave function of the condensed boson, (pc is defined as a solution of 
the following Gross-Pitaevskii equation: 



(16) 



2 

where ph is the boson chemical potential. In the ground state the fermion phase-space distribu- 
tion function is given by the Thomas-Fermi (TF) approximation 

f{r,p) = e[iif-e{r,p)l (17) 

with 

e{r,p) = -^p' + UF{r), (18) 

where is the fermion chemical potential. In this TF approximation the fermion density /9_p is 
obtained as the solution of the following equation: 

^ {67rVF(r-)}'/' + \msuo){rl + ^z^) ^ hBFPB{r) = Pf. (19) 

Here we iterate solving the boson wave function with Eq. (jl6p and searching the fermi energy fj,f 
in Eq. ()19p to give the correct fermion number. 
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In Fig. [T]we show the density distribution of the boson and fermi gases with the boson-fermion 
couphng hBF = —QBE hsp = Qbb (b) and hsF = "^Obb (c). The sohd and dashed hnes 
represent the results of boson and fermi gases, respectively. The density distribution of the 
bose gas is negligibly changed while the fermion density distribution in the boson occupation 
region varies as the boson-fermion coupling increases. The fermion density distribution is central 
peaked when Hbf = —gBB, hat when hsF = Sbb, and surface peaked when hsF = "^Obb- 
This boson-fermion coupling dependence of the fermion density distribution can be easily ex- 
plained as follows. When the boson number, N^, is very large, the boson density distribution 
can also be given by the TF approximation as 



PB{r) ^ 



- \{fT + - hBFPF{r) 



(20) 

QBB V -l - ^ " ' ■ ■ 'J 

Note that the boson density in Eq. (l20p is defined within a region, and it is zero outside of that 
region. 

Substituting Eq. ([20|) into Eq. (fT9]) . we can get 

1/2 ^2/3 h\p Kbf 1 / 2 ^Bf\ ^2 /oi x 

(67r pf) PF = P-F PB - - rufUJf ^ , 21) 

2m/ gBB gBB 2 V gBB J 

where = r'^+n'j^z'^. From this equation we can easily know that pp = const, when hsF / gBB = 
1 and pB > 0. Furthermore the derivative of the fermion density, pp, with respect to is given 
as 

The stability condition 



18t |. that {dph/dpB){dpf /dpp) — {dpb/ dpF){dpf / Bpb) > 0, restricts 
the value of the fermion density as pp < 47r^(7^g/3mj/i^p . In the boson occupation region, 
therefore, the derivative of the fermion density is negative, dpp/d^ < 0, when hBp/gBB > rnfLOp 
and positive, dpp/d^ > 0, when hBp/gBB > nifuj'j. 

As mentioned before, we take the parameters to be mfUi'j = 1 and = 1 in the present 
calculation. Thus the TF approximation can explain the relation between the boson-fermion 
coupling and the fermion density distributions. 



B. Dipole Oscillation 

In this subsection we show our actual results of the numerical simulations on the dipole os- 
cillations. Here we define the center of mass (CM) position on z-coordinates for bosons and 
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fermions as zb and zp, respectively. We will discuss the oscillation behavior by examining the 
time-dependence of and zp- 

In actual simulations we boost the condensed boson wave function and the fermion test-particles 
at the starting time, r = 0, in the following way: 

<Ac(r,r = 0) = e*"^^#(r), (23) 

Pz{i)=pi^H^) + mfUjfXF (24) 

with the boost parameters and Xp, where the superscript {g) represents the wave function 
and the coordinates of the ground state. These transformations give the current density of 
boson, jsif), and fermion, jp{r), as 

jB{r,r = 0) = p'^i\r)XBZ, (25) 
j^(r,r = 0) = pf{r)u;fXpz. (26) 

If hBF = 0, the time dependences of zb,f become zb = ABsin(T) and zp = Xp sm{ujfT); the 
boost parameters and Xp correspond to the initial amplitudes of zb and zp, respectively. 
In Fig. [2] we show the time-dependence of the zb,p when the boson- fermion coupling Ubf = 
—gBB (a), hBF = —0-5gBB (b), hBF = Qbb (c) and Hbf = 2gBB (d). The dashed and solid lines 
represent the results of zb and zp, respectively. In all calculations we choose the out-of-phase 
at the beginning between the boson and fermion oscillations by taking the initial condition to 
be Xb = 0.4 and Xp = -0.4. 

While the boson oscillations are monotonous, the fermion oscillations have damping in the early 
time stage. When Hbf = ^qbb, the amplitude of the fermion oscillation is about 0.4 at the 
beginning, but it decreases and becomes about 0.1 — 0.2 after the damping (r > 60). This 
damping becomes slower when the coupling is weaker, Hbf = —0-5gBB (Fig-Eb) • Furthermore, 
we see that after the damping zb and zp oscillate with the same period, and their relative phase 
becomes in-phase when Hbf < and out-of-phase when Hbf > 0. Note that, when Hbf = '^9bb, 
the damping of zp is not clearly seen; the damping is too fast, and the period of the oscillation 
becomes the same with that of zb in the early time stage. 

In order to confirm the above comment about the relative phase after the damping, we calculate 
the dipole oscillation with Hbf = Qbb, which is started with in-phase at the beginning (A^ = 
Xp = 0.4). The results are shown in Fig. [3l The relative phase between the boson and fermion 
oscillations becomes out-of-phase in the later time stage after the damping. In order to clarify 
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it, we plot the same quantities shown in Fig.[2h, Fig. and Fig.[3]in later time, 120 < r < 150, 
in Fig. m We can see that zp oscillates almost with the same period of zb- 
As seen in Figs. [2]-|31 the behaviors of the fermion oscillations are not so simple and imply that 
the fermion oscillation includes various modes. In order to inspect these properties more, we 
calculate the strength function defined as the Fourier transformation of zb^f'- 



oscillations with Ubf = —gsB for the bosons (a) and the fermions (b) and with Ubf = —0-5gBB 
for the bosons (c) and the fermions (d) in Fig. [5l and those with Hbf = Qbb for the bosons (a) 
and the fermions (b) and Hbf = '^dBB for the bosons (c) and the fermions (d) in Fig. [H 
First we note that the boson strength functions have only one sharp peak, which is consistent 
with the monotonous behavior of the boson oscillations. Hence the frequency of the peak position 
can be considered to be the intrinsic frequency of the boson oscillation; we define u;^ as this 
frequency. Second the fermion strength has three peaks. One peak appears at a; = and 
another peak appears at the trapped frequency, lo = ujf ^ 1.47 in all cases, whereas the position 
of the other peak depends on the coupling constant. These three peaks must correspond to 
certain modes. For convenience we refer to these three modes as mode-1, mode-2 and mode-3 
in order, and define uj^ as the frequency of the mode-3. 

The sign of the fermion strength function at w = uj^ (mode-1) is plus when Hbf > and minus 
when Hbf < 0, while the signs of the strengths at u! = ujf (mode-2) and u! = uj-^ (mode-3) are 
minus in all oscillations. In Fig. [71 furthermore, we plot the strength functions of the bosons (a) 
and the fermions (b) with Hbf = Bbb using the in-phase initial condition (Ab = Xf = 0.4). We 
see that the signs of the strength functions at w = and uj = become plus. 
In this initial condition the sign of the boson strength function is always plus. The sign of the 
strength exhibits the relative phase between each mode and the boson oscillation; the plus and 
minus signs show the in-phase and out-of-phase, respectively. Then the above results imply that 
a choice of the phase at the beginning determines the phases of the mode-2 and the mode-3, 
while the phase of mode-1 is determined by the boson-fermion coupling. These results exhibit 
typical behaviors of the forced vibration with the external force caused by the boson oscillation, 
which is discussed in the following subsection. 




(27) 



In this work we fix that tj = and tf = 200 {Q^ ). We show the strength functions of the dipole 



9 



C. Modes of the fermion oscillation 

In this subsection we examine the above three modes for the fermion oscillations by performing 
the test simulations shown in the following. 

First we calculate the fermion dipole oscillations with the boson motion frozen; namely the 
fermions move in the fixed potential Up^r) in the ground state. In Fig. [8] we show the results 
with the boson-fermion coupling hsF = —QBE (a), hsF = 933 (b) and hsF = "^933 (c). In each 
result we see a beat and damping: the amplitude of zp oscillation is about 0.4 in the beginning, 
and it becomes about 0.05 after the damping. 

In Fig. [9] we plot the strength functions of these fermion oscillations with the boson-fermion 
coupling hsF = —933 (&)> hsF = 933 (b) and hsF = '^933 is). In all the results there are two 
peaks ai oj = ujf and uj = oj^^, but no peak at w = w^. 

These results demonstrate that mode-2 and mode-3 are intrinsic modes for the fermion oscilla- 
tion, and that the mode-1 is caused by the boson oscillation. A mixture of mode-2 and mode-3 
arises the beat, and this mixing and their widths make the damping in the fermion oscillations. 
Next we simulate the dipole oscillation with the initial condition Xb = 0.4 and Xf = 0. In Fig.fTOl 
we show the time-dependences of zb and zp with the boson-fermion coupling Hbf = —933 (a)) 
hBF = 933 (b) and Hbf = "^933 (c). In Fig. [TTl furthermore, we plot the strength functions of 
the boson oscillation with Hbf = —933 (a) and with Hbf = 933 (b) and those of the fermion 
oscillation with Hbf = —933 (c) and with Hbf = 933 (d). In the fermion strength functions 
(c,d) there are two peaks at u> = uj^ (mode-1) and u> = uj^ (mode-3), and no peak at uj = ujf 
(mode-2). In this initial condition, Xp = 0, the boson oscillation is a trigger of the fermion 
oscillation and moves only fermions inside of the boson occupation region. Hence these results 
demonstrate that mode-2 and mode-3 are the dipole motions contributed from the fermions 
outside and inside of the boson occupation region, respectively. 

The condensed boson density distributes in a smaller region than the fermion one, so that 
the fermion potential. Up, is separated into two regions, the inside and outside of the boson 
occupation region. Outside the boson occupation region the fermion potential. Up, has a simple 
harmonic oscillator shape with the trapped frequency coj. Inside the boson occupation region, 
furthermore, the fermion potential around the central region can also be approximately described 
as a harmonic oscillator potential with a different trapped frequency. 

Thus this fermion motion can be described with two kinds of fluids corresponding to the above 
two modes. One fluid moves outside the boson occupation region, and makes a dipole oscillation 
mode with the frequency, to = coj. The other fluid moving in the inside region makes the 
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dipole oscillation with uj = . However the fermion potential does not have a simple harmonic 
oscillator shape in the boson occupation region, and then their dipole motions, particularly for 
the latter fluid, are not harmonic, and the oscillation amplitudes decrease. 
Now we examine the above consideration in quantum calculations. We calculate the excited 
states in RPA and compare the results with those in our time-dependent approach. For this 
purpose we get the fermion wave functions in the ground states in the Hartree-Fock (HF) approx- 



281]. In the HF calculations the fermion 



imation and solve the RPA equation in the way of Ref. 
number, A'^^, must be dependent on the subshell closure of fermion single particle states, and 
changed from Nf = 1000; for example Nf = 989 with hsF = —Obb, but this slight change 
hardly affects the final results. We use seventeen particle-hole states for bosons and about five 
hundred ones for fermions. 

Then we calculate the transition amplitudes from an excited state | > with excitation energy 
iVn to the ground state |<I>o > for bosons, Ab, and for fermions, Ap, which are written as 

AsiiOn) = < ^-nl^Bl^O >, (28) 

ApiUn) = <<^n\ZF\<^0> (29) 

with 

= j d^r<p*{r)z4>{r), (30) 

Zf = [ Sri)*{r)z^{r). (31) 



In Fig. [12] we show the boson Tb(u;„) = |ylB('^n)P/4'7r (a) and fermion transition strengths 
Tp{uJn) = |^F(<^n)P/4vr (b) as functions of the excitation energy uj with Ubf = —9bb- 
In order to visualize the distribution for the fermion transition strength, we introduce an artificial 
width, r, and draw the curve of Tp{uj) defined by 

r|H = ^E^H^„) p^(J_,„). . (32) 



n 



In our calculation we take F = 0.01, and plot 0.018 x Tp{uj) in the second panel (b), where the 
arbitrary factor, 0.018, is introduced to plot the curve with the same scale of Tp in the same 
figure. 

In these results there is a clear peak at cj = 1.0 in the boson transition and two clear peaks 
at w = 1.50 and 1.61 in the fermion transition. These two peak positions of Tp are close to 
the frequency of the mode-2, w/, and that of the mode-3, a;{, in our time-dependent approach. 
Furthermore we see several modes with the excitation energy between ojf and w;^. 
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We would like to comment that there is a peak Tp = 0.075 at liJ ~ 1.0, but it is too small to be 
visible. In the RPA calculation, thus, there is not any any strong co-moving mode, and hence 
the fermion transitions must be almost explained as single particle processes. In Fig. [T2b we 
plot the fermion transition strength, Tp, in the single particle process; the results of Figs. 112b 
and [T2b are almost the same. 

Hence the dipole oscillation behavior can be explained as follows. The condensed bosons oc- 
cupy one single particle state, and the boson oscillation has one mode and does not damp. In 
contrast the fermions occupy many single particle states, and the fermion oscillation includes 
two intrinsic modes which are contributed from the fermi gases outside and inside of the boson 
occupation region, respectively. Since the fermion potential cannot be critically separated, there 
are several modes with frequencies between ujf and w^. Because of this statistical difference, 
furthermore, the boson oscillation almost one-sidedly affects the fermion oscillation. Then the 
fermion oscillation has one more mode caused by the external vibration force which is the boson 
oscillation. As a result the fermi gas gradually loses its intrinsic modes and finally oscillates 
with the same frequency of the condensed bose gas. 



D. Comparison with the Sum-Rule Approach 

In this subsection we discuss the boson- fermion couplin g de pendence of the oscillation frequencies 
by comparing our results with the sum-rule approach 29(1 . 

In Fig. [11] we show the intrinsic frequencies of boson and fermion with full diamonds and circles, 
respectively. For references the fermion frequencies with the boson motion frozen are plotted 
with the open circles which are connected with the long-dashed line, and the trapped frequency 
of the fermion with the dotted line. In addition we give the results of the RPA, which is defined 
as the peak energy of Tp, with the asterisk. Our results about the fermion oscillation in the full 
calculation (full circles) agree with those in the boson motion frozen (open circles) and those in 
RPA (asterisk). 

□ 

The sum-rule approach [29] gives the intrinsic frequencies of the boson and fermion dipole 
oscillations as 



= j^^KP, (33) 
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with 

VnP = hsJSr'-^'-^. (35) 
J oz oz 

In the same figure we also plot the results of the sum-rule approach with the thick dashed lines. 
The above results in the sum-rule approach agree with our results when hsF < 0, though the 
two results show discrepancy when hsF > 0. 

It is known from Eqs. (j33p . (j34p and (j35p that the frequency calculated with the sum-rule is sen- 
sitive to the density distribution. As discussed in Sec.III-A the derivative of the fermion density 
with respect to the radial coordinate in the boson occupation region is negative, dpp/dr < 0, 
when hsF ^ Obb, and positive, dpp/dr > 0, when hsF ^ gsB, while the derivative of the 
boson density, dps/dr < 0, in the both coupling region. As the boson-fermion coupling, hsF, 
increases, hence, the fermion frequency in the sum-rule approach decreases when Hbf < 0, is 
not varied when < /ib_f < 1, and increases when hsF > . 

The sum-rule approach is available only when the transition strength is concentrated to one 
excited state. In BF mixtures, however, the fermion transition strength distributes mainly to 
the two modes, and then the sum-rule approach cannot predict the correct intrinsic frequency 
of the fermion oscillations. 

Both our time-dependent approach and the RPA calculation explain the two modes in the 
fermion oscillation. However the RPA can be available only in minimal oscillations, and we 
should examine the time dependence of zb,f in RPA. 

Here we calculate the time evolutions zb and zf in RPA. By using the initial boost given in 
Eqs. (l23p and (j24l) with the perturbative way, the initial state in the dipole oscillation becomes 

|$(r = 0) > = eH-^fl^i'+-^f'^/}|$o > 

^ |$o > +i {Afiifc + Af^/} |^>o > • 

= |$o >+iYj {Ab^bK) + Af^fK)} |$n > • (36) 

n 

This initial boost is equivalent to that given in Eqs. ()23p and ()24p in the first order of Xb,f- 
Then, the time-dependence oi zb and zf are given as 

n 

= —'^AB{uJn){XBAB{uJn)+XFAF{uJn))sm{u;nT), (37) 

n 

f n 
2 V - 

= — > ^ AF(w„)(AB^B(a;„) + AF^F(t^n))sin(w„r). (38) 

iV f ^ — ' 

n 
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In Fig. [T3]we show the results for boson (dotted line) and fermions (solid line) in the out-of-phase 
oscillations with hsF = —QBE (a) and hsF = dBB (b) by using the condition Xb = —Xf = Az. 
The amplitudes of the oscillations are scaled by Az. Comparing these results with those in the 
time-dependent approach (Fig. [2]), we can see a clearer beat and weaker damping in the fermion 
oscillations, and the period of the fermion oscillation dose not become the same as that of the 
boson oscillation even in the later time stage. The boson oscillation does not affect the fermion 
oscillation in RPA as strongly as in the time dependent approach. 

In order to inspect the reason for this difference more, we show the density distribution and the 
velocity fields at r = (a) and at r = 18.5 (b) in Fig. [T5j The dashed and solid contour lines 
represent the density distributions for boson and fermion, respectively, and the arrows indicate 
the fermion velocities. At both times the z-component of the CM position of fermi gas is zero, 

ZF = 0. 

At r = the boson and fermi gases move with unique velocities (a) in all positions. At r = 18.5, 
however, the velocity of the fermi gas is directed along the contour lines of the density, and its 
density becomes dilute in the space region, z < and {rxl ^ 1. The fermi gas moves to the 
surface part of the bose gas because of the repulsive force from the bose gas. This motion is 
very natural from the aspect of the hydrodynamics. 

The RPA approach can describe only minimal oscillations with one-particle one-hole excitation 
energy, where the density distribution is slightly changed from that at the ground state. In the 
case of Hbf = Obb and the boson motion frozen, for example, its excitation energy is about 80 
KIb, and the system cannot be described with simple one-particle and one-hole states. In such 
non-small amplitude oscillations the fermi gas changes its density distribution and moves to the 
place with low resistance. Hence the actual oscillation of the BF mixing gases are not so simple 
to describe with the RPA and sum-rule approaches. 

IV. SUMMARY 

In this paper we study the collective dipole oscillation in the BF mixtures by solving the time 
evolutions of the system directly with the TDGP and Vlasov equations. The calculational 
results reveal that the boson and fermion oscillations make different behaviors reflected by their 
statistics. The condensed bosons occupy one single particle state, but the fermions distribute 
in many single particle states. Then the boson oscillation has only one mode and does not have 
damping, but the fermion collective oscillation includes various modes. In the BF mixtures, 
especially, the fermi gas is separated into the two regions, inside and outside of the boson 
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occupation region, and these two fluids oscillate with diff'erent periods. 

Furthermore the boson oscillation almost one-sidedly affects the fermion oscillation, and then 
the fermion oscillation includes three modes: one is a mode of the forced vibration caused by 
the boson oscillation, the other modes are two fermion intrinsic modes contributed from the 
fermi gases inside and outside of the boson occupation region, respectively. In addition there 
are several small modes in the fermion oscillations with frequencies between the two intrinsic 
frequencies. These modes cause the beat and damping in the fermion oscillation. As the result 
the fermion motion gradually loses the strength of its intrinsic modes and finally oscillates with 
the same frequency of the condensed bose gas. 

The RPA approach can also explain the frequencies of the boson and fermion intrinsic modes. 
The qualitative behavior of the transition strength in RPA oscillations is similar to the strength 
function in our approach. Nevertheless the time evolutions of zp are not the same, particularly 
in the later time stage, because the forced oscillation modes, mode-1, does not largely contribute 
to the whole oscillation in RPA. RPA is available only in the minimal oscillation, and cannot 
describe the change of the density distribution in the time evolution process. In actual exper- 
iments the amplitude is not so small, and hence we must solve the time-dependent process of 
the collective oscillations directly. 

In this work we do not take into account two body collisions and thermal boson 
system ^ Nf at zero temperature, the thermal bosons are very few, and then the two 
body collisions do not play significant roles in the dynamical process. However, it is not easy 
to establish gases at such low temperatures in actual experiments. In the future we need to 
introduce the two body collision terms into our approach 
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FIG. 1: The density distribution of the BF mixing gas with the boson- fermion couphng h 
{&), hsF = 9bb (b), and Hbf = gsB (c). The dashed and sohd hnes represent the density ( 
and fermi gases, respectively. 
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FIG. 2: Time evolution of the center-of-mass position of the boson (dotted Hnes) and fermion (soHd Hnes) 
with the boson-fermion coupHng hsF = —gsB (a), hsF = —O-^gsB (b), hsF = gsB (c) and hsF = '^gsB 
(d). The initial condition is taken to be As = 0.4 and Xf = —0.4. 
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FIG. 4: Time evolution of Zb (dashed lines) and zp (sold lines) with the boson-fermion coupling Hbf = 
—Qbb (a), hsF = Qbb (b) and Hbf = Qbb {c). In the first and second panels the oscillation is started 
with the out-of-phase at the beginning, and the last ones with the in-phase. 
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FIG. 5: Strength functions of the boson oscihation in upper panels (a,c) and fermion oscillation in lower 
panels (b,d) with Hef — —gsB (a,b) and with Hef = —O-^gsB (c,d). 
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FIG. 6: Strength functions of the boson oscillation in upper panels (a,c) and fermion oscillation in lower 
panels (b,d) with Hbf =933 (a,b) and with Hbf = ^933 (c,d). 
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FIG. 7: The strength fmiction of the fermion dipole oscillation with hsF = Qbb, but using the in-phase 
initial condition (As = Ai? = 0.4). 
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FIG. 9: Strength functions of the fermion oscillation with hsF = 9bb (a), 0-5gBB (b), —0.5gBB (c) and 
with Hbf = —Qbb (d) when the boson motion is frozen. The dotted lines indicate the fermion trapped 
frequency ojf. 
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FIG. 10: Time evolution of zb (dotted lines) and zp (sold lines) without the boson boost at the beginning 
with the boson-fermion coupling hsF = —Qbb (a), ftsF = Qbb (b) and hBF = "^Qbb (c). 
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FIG. 11: Strength functions of the boson (a,c) and fermion (c,d) osciUations with hsp = —gsB (a,b) and 
with hsF = 9bb (c,d). The details are shown in Fig. [TOl 
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FIG. 12: Transition strengths versus tlie exeitation energy for tlie boson (a) and ferniion (b) oscillations 
with hsF = ^9bb in RPA. The fermion transition strength in the single particle process are also ploted 
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FIG. 13: Intrinsic frequencies of boson and fermion oscillations versus the boson-fermion coupling. Full 
diamonds and full circles represent the intrinsic frequencies of boson and fermion oscillations, respectively. 
The open circles, which are connected with the long-dashed line, indicate the fermion intrinsic frequency 
with the boson motion frozen. The thick dashed and dotted lines denote the results of the sum-rule and 
the trapped frequency of the fermion, respectively. 
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FIG. 14: Time evolutions of zb (solid line) and zp (dashed line) in RPA with hsp/gBB = —1 (a) and 
h-Bp/gBB = 1 (b). 




FIG. 15: Density distribution and velocity fields, at r = (a) and r = 18.5 (b). The dashed and solid 
contour lines represent the density distributions of the bose and fermi gases, respectively, and the arrows 
indicate the velocities of the fermi gas. 



